Point-Contact Conductances from Density Correlations 
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We formulate and prove an exact relation which expresses the moments of the two-point conductance 
for an open disordered electron system in terms of certain density correlators of the corresponding 
closed system. As an application of the relation, we demonstrate that the typical two-point con- 
ductance for the Chalker-Coddington model at criticality transforms like a two-point function in 
conformal field theory. 
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In mesoscopic physics, as well as in quantum physics 
at large, open systems are distinguished from closed ones. 
The difference is most evident from the nature of the cor- 
responding energy spectra: when a closed and finite sys- 
tem is opened up, the discrete set of stationary solutions 
of the Schrodinger equation turns into a continuum. In 
the open case, one usually looks for solutions in the form 
of scattering states, where the wave amplitudes are pre- 
scribed in all incoming channels and the outgoing waves 
are then related to the incoming ones by a linear operator 
called the scattering matrix or ^-matrix. 

One may ask whether one can connect the scatter- 
ing matrix of an open system with the properties of its 
closed analog. In spite of the physical distinction between 
open and closed systems, such connections do exist. For 
example, for nonintcracting quasi-one-dimensional elec- 
trons subject to weak disorder it has been found Q that 
the parametric correlations between the eigenphases of 
the ^-matrix exhibit behavior similar to that of the en- 
ergy levels of the closed system. For another example, 
the jS-matrix (or rather its extension including evanes- 
cent modes) has been utilized pj in numerical algorithms 
that compute large numbers of energy eigenvalues for a 
chaotic closed billiard. An intimately related scheme is 
provided by Bogomolny's semiclassical T-operator j|. 

All these theoretical developments were primarily con- 
cerned with the energy levels and their statistics. In the 
present Letter we address a different question: suppose 
a closed system has been opened, what relations, if any, 
exist between the statistical properties of the S-matrix 
and the stationary wavefunctions of the closed system? 
This question is of practical relevance, as the S-matrix 
determines the electrical conductance of an open meso- 
scopic device. It will be seen that, provided a few well- 
formulated conditions are satisfied, an informative an- 
swer can be given. While this seems unlikely if the sys- 
tem is perturbed strongly, say by making the boundaries 
transparent, the situation is different if one uses point 
contacts which establish the connection to the outside 
through ideal one-dimensional channels. For concrete- 
ness and simplicity, we are going to develop the theory 
for the case where time is a discrete variable, but we 
anticipate that a variant of the relation to be derived re- 



mains valid for systems with continuous-time dynamics 
generated by a Hamiltonian. 

We begin by briefly reviewing the concept of network 
models. Wave or quantum-particle propagation in a ran- 
dom medium or a chaotic cavity can be modeled by a 
network (or graph) consisting of internal channels repre- 
sented by links and of local scattering centers represented 
by nodes. For a network with N links, the state is speci- 
fied by N complex amplitudes that assemble into a vector 
$ = {ipi}i=i,—,N- The network evolves in time by re- 
peated application of a unitary matrix U G \J(N), which 
propagates states between integer times t: ty(t + 1) = 
U^(t). Stationary states tyj = {ip^i} and quasi-energy 
levels < lui < . . . < < 2n are constructed by solv- 
ing the eigenvalue problem U^j = e^&j. By the local 
intensity pjj of we mean the squared amplitude l^,;) 2 
of the normalized state $ j at link I. 

The previous definitions apply for closed networks, 
where each link starts and ends at a network node. To 
define a conductance, however, one needs to connect the 
network to external charge reservoirs. We do this by in- 
troducing two interior point contacts, i.e. by cutting the 
network open at two links / and m, thereby producing 
two incoming and outgoing channels. If we imagine these 
channels to be connected to Ohmic contacts, the upshot 
of the surgery done on the network is to define a two- 
terminal conductance between the points I and m. We 
denote it by T[ m and refer to it as the two-point conduc- 
tance for short. The inset of Fig. [I] shows the basic setup 
at the example of the Chalker-Coddington network 0] . 

Given the 2x2 scattering matrix S = ( g lm j 

which consists of the amplitudes for transmission and re- 
flection by the open network, the Landauer-Biittiker for- 
mula states that Tj m (at zero temperature, and in units of 
the conductance quantum e 2 /h) equals the transmission 
probability |Si m | 2 . For the following it will be crucial 
that 5*, though primarily defined for the open network, 
can also be expressed || entirely in terms of the unitary 
time-evolution operator U of the closed system: 

S aP = (a\{l-UPiP m )- x U\l3) (a,/3 = l,m) (1) 

where = 1— Expansion of (1 — UPiPm)^ 1 in a 

geometric series shows that S a p is a sum over path am- 
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plitudes for paths from (3 to a. The projection operators 
Pi and P m serve to truncate paths as soon as they arrive 
at the outgoing channel of either terminal I or m. 

Via relation (|]), every disorder average (• • •) of the 
closed system determines a disorder average pertaining 
to the S'-matrix of the open system. In particular, the av- 
erage of any function of the two-point conductance is de- 
fined, by (F(T lm )) = (f (|(Z|(1 - U P l P m )- 1 U\m)\ 2 ')) . 

Starting from some disordered network model, our aim 
in the sequel is to relate the distribution of two-point con- 
ductances to the statistics of stationary wavefunctions. 
As will be seen, the relation we are going to derive does 
not depend on the details of the model chosen. The only 
requirements we need are two conditions constraining the 
disorder average: the phases of the reflection coefficients 
Su and S mm must be (i) uniformly distributed and inde- 
pendent of each other, and (ii) statistically independent 
of Tim . These conditions are automatically satisfied for 
the Chalker-Coddington model and more generally for all 
models with local phase invariance. 

Given some choice of model with disorder average (•••), 
we introduce an additional average (• • ■) for observables 
that can be expressed by the stationary states VE^-: 



(A) c := - I du;^ - "i)M*j)) > (2) 
-"Jo j=l 
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where v is the level density at u> = 0, and e sets the size 
of an energy window. In what follows, we take e — > 0, 
restricting the sum over states to the first level, io\. 

We are now in a position to propose a statistical re- 
lation between Tj m and the intensities pi = IV'i.zl 2 an d 



l^i 



Let f(x) be any function defined for x in 



the positive real numbers. Then we claim Q 

2tw (p m f {p m /pi)) c = (F(Ti m )) , where (3) 
F(T) := J* g / (T- 1 !! - c^VT^T?) ■ (4) 



A first remark on (||) is that, in view of the scaling v ~ N 
and p m ~ A -1 (on average over the statistical ensemble), 
the system size L d ~ A drops out of the combination 
v Pmf{Pm/ Pi), reflecting the fact fjj that Ti m tends to a 
finite limit as L is sent to infinity. Note that the plain 
correlator v 2 {pip rn ) does not share || this property of 
size independence, which is why any attempt to express 
it through such relations as (|[^) is met with failure. 

More explicit relations are obtained by adopting spe- 
cial choices for /. Of particular interest is an expression 
giving the typical conductance T* m p = exp (lnTi m ): 



(lnT im ) = 2-kv (p m ln(pi/p m )) e , 



(5) 



which follows from (J3j|4|) on taking f(x) = — lnx. An- 
other simple relation results on setting f(x) = (1 + x)~ 1 : 



(VTm) = AlTV (pipm/(pl + Prn)) c ■ 

Furthermore, by inverting the relation gotten by setting 
f{x) = x n , and analytically continuing from integer n to 
complex z, it can be shown with some effort that 
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T(k + 2 - z)T(l - k - z) 
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{Pm Pi ) c 



Note that the sum over k is terminated by the poles of 
the gamma function when z is a negative integer or zero. 
For all other values of z, the series will in general be 
divergent, since the polynomial decay of the product of 
inverse gamma functions does not suffice to control the 
growth with k of the moments (p^ 1 p^ k ) . Fortunately, 
if z = n + 1/2 with n a positive integer, the series can be 
rearranged to produce a finite sum involving polynomials 
in the bounded random variable pip m /{pi + Pm) 2 < 1/4: 
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(6) 



where p n = 2 An + 2 vv({p lP m) n+1 {pi + p™)- 2 "- 1 )^ We 
will put this formula to good use below. 

Our main result is relation (^). Its proof is based on 
the observation that for a network at resonance, i.e. for 
u>i = 0, the corresponding eigenvector of U can be 
interpreted as a scattering state of the network opened 
at the links / and m. In other words, the two amplitudes 
?/>i,z and t/)i t m of an eigenstate ^i = U^i constitute at 
the same time an eigenvector ("0i,z, V'l.m) with eigenvalue 
unity of the scattering matrix 5*. Thus we have Sui^i^i + 
Si m ipi,m = ^1,1 and therefore, assuming ^i,jS; m ^ 0, 
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(7) 



The next step will be to take the disorder average of 
this equation, viewing both sides as functions on the 
closed network. Of course, since (0) holds only for net- 
works at resonance, it is imperative that we restrict the 
disorder average to the resonant realizations. We do this 
by tuning the random phase tp m at link m to a function 
(p m = ifimiO °f au other random variables £ = {£i, ^2) ...} 
so that <Pm(Q) = is identically satisfied for all £ 
H . The restricted or resonant disorder average thus de- 
fined is denoted by (• • -) rcs m . By applying some function 
/ to the identity (Q) and averaging in this manner we get 
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where we have set (pi = arg Su. Now recall that, accord- 
ing to conditions (i) and (ii), the phase (pi is statistically 
independent of ip m and T; m and is uniformly distributed. 
Consequently, the preceding equation continues to hold if 
we integrate the r.h.s. with measure dipi/2-K over [0,27r], 
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thereby turning it into (F(Ti m )) TCS m with F defined by 
(||). Moreover, condition (ii) entails that the resonant av- 
erage (F(Ti m )) rcs m actually equals the natural average 
(F(T lm )). Thus 



(f(Pm/Pl)) Ies ,m = ( F i T lm)) ■ 



(8) 



Next we transform the average over resonant networks 
on the l.h.s. to the eigenfunction average (||), as follows: 

2icv (A) c = hm — / duj - w)A(^ 1 )) = 
£ Jo 



2,T 



d<p m (£(wi(£,¥>m))^(¥i))< 



A(gi) 
4 \|e?wi/ctyj m | 



Setting A 



f(p m /pi) and using (rah we arrive at 



(F(T lm )) = 2-ki 



duo i 



d(p„ 



f 



Pi 



Finally, since a phase variation 5ip m perturbs the time- 
evolution operator U by i(S<p m )\m) (m\U , standard first- 
order perturbation theory shows that the level velocity 
\dwi/d(p m \ equals p m . This concludes the proof. 

It is easily seen that the above line of reasoning does 
not depend in any essential way on assuming a network 
model with discrete time. In fact, a slight variant of 
the argument goes through for continuous-time dynam- 
ics generated by a Hamiltonian H . What one needs to do 
is attach two identical and strictly one-dimensional leads 
I and to to the system, which are then either connected 
to charge reservoirs, or else closed off by imposing reflect- 
ing boundary conditions at variable lengths L; and L m . 
The intensity pu m ) 01 the closed system is interpreted as 
j*o+ ^ x |^,^|2 w h ere \ is the wavelength and Xq, xo + X 
are in the closed lead / (to). The normalization factor 
in Eq. (||) has to be replaced by (p m ) _1 . To tune the sys- 
tem to a given resonance energy E, we adjust L m . 

At the moment it is not clear to us to which extent, 
if at all, the relation applies to systems with more than 
two point contacts or more general types of contact. This 
point deserves further investigation. 

When states are localized, relation (||) implies that 
the decay length of the typical two-point conductance 
agrees with the typical localization length of the states 
of the closed system. While this has long been known 
for strictly one-dimensional systems po[ , it is now seen 
to hold exactly under the much more general conditions 
we have specified. By word of caution, however, we recall 
that our formula does not cover correlators like v 2 (pip m ), 
so it does not settle a recent debate as to whether two- 
scale localization in a magnetic field can occur for 
closed systems when it does not ]T^ ] for open ones. 

We now demonstrate the utility of the result (S) by 
applying it to a challenging problem: critical transport 




FIG. 1. Half-integer moments (T' I+1/2 ) for n = 0, . . . , 29 
of the two-point conductance of the Chalker-Coddington 
model at criticality. The four data sets correspond to con- 
tact distances x = 5, 10, 15 and 20 (top to bottom) in a 
cylinder geometry with length L — 100 and circumference 
W = 10. Solid lines represent the theoretical prediction. 
Inset: Chalker-Coddington network connected to external 
charge reservoirs via two interior point contacts. 

at the integer quantum Hall (IQH) transition p3[ . To 
make analytical progress with that transition, one needs 
to confront the question what can be learned from con- 
formal field theory. For critical systems without disor- 
der, the principle of conformal invariance has proven to 
be very powerful, by posing strong constraints on cor- 
relation functions and determining how they transform 
between various geometries [{lif . The question whether 
the same principle applies to 2d disordered critical points 
was first raised by Chalker [fisfl . The answer turned out 
to be negative for the correlation functions he considered. 
However, as we now proceed to demonstrate, conformal 
invariance does govern the two-point conductances. 

An analytical formula for the distribution of two-point 
conductances T = XJ m for critical quantum Hall systems 
in the infinite plane was proposed in ||. Given that 
formula, conformal field theory predicts the probability 
measure d/i(T) for the conductance between two points 
with coordinates (0,0) and (x,y) on an infinitely long 
cylinder of circumference W to be given by 



dfi(T) = piT- 1 - l)T- 2 dT 

27r -l/2 c -X t /4 



(A t lnC) 3 / 2 



In Qftft 

y sinh 2 t — p 



(9) 
(10) 



arsinh^/p 

£ = (W/na) |sinh(7r(x + iy)/W)\ , 

where X t is a critical exponent (unknown as yet from 
theory), and a is a nonuniversal microscopic scale that 
sets the length unit in which to express the correlation 
functions of the particular IQH model considered. 

To test this prediction, we numerically calculated the 
eigenfunctions of the critical time-evolution operator U 
of the Chalker-Coddington model, in a cylindrical geom- 
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etry with large aspect ratios L/W, 100 < L < 200 and 
3 < W < 30. An ensemble of 5000 disorder realizations 
was employed. Using Eq. ([|) we computed the moments 
(T ,l+1 / 2 ) by averaging over the pairs of links (l,m) with 
a given distance \x\ and their difference vector pointing 
along the axis of the cylinder. The great advantage of this 
indirect approach (in contradistinction with the direct 
method of S) is that a large number ~ N 2 of data points 
can be extracted from a single eigenfunction 'J. Since the 
numerical effort of calculating W and T; TO is comparable, 
we gain in efficiency by a factor iV 2 ~ (LW) 2 , which can 
be many orders of magnitude. 

In Fig. [l] theory and numerics are compared for W = 
10, L — 100, and \x\ = 5, 10, 15, 20 (in plaquette units). 
The theoretical curves were obtained by numerically in- 
tegrating (||) against T n + 1 / 2 , We find that an exponent 
X t = 0.54 ± 0.01 (and a = 0.20) yields the best fit. The 
statistical errors are estimated to be comparable with the 
symbol size. The agreement is clearly very good, showing 
that the result predicted by conformal field theory inter- 
polates correctly between 2d behavior at short distances 
and quasi- Id behavior at large scales. 

The most clear-cut test of conformal invariance is of- 
fered by the typical two-point conductance. By integrat- 
ing (||) against InT, one finds a very simple law j^]: 

T £o);W = 10*7™) sinh(nx/W)\- X * . (11) 

Fig. ^ plots (InT) versus ln|x| on cylinders of length 
L = 200 and circumferences W ranging over one decade, 
from 3 to 30. We here used 2800 disorder configura- 
tions for each W. In the inset (InT) + XtlnVF is plot- 
ted against In |x/W| with an exponent X t = 0.57, deter- 
mined with an error of 0.05 by optimizing the match of 
the numerical data with the prediction ( pi] ) (a = 0.20). 
It is seen that, within the statistical errors, the data 
collapse onto a single curve given by (|il|). This means 
in particular that the short-distance algebraic behavior 
T typ = \x/a\~ Xt crosses over to the exponential law 
T typ ex e _A l x l for |x| ^> W, with the Lyapunov expo- 
nent being A = irX t /W. While the relation XW = const, 
has been checked before fl6)| , the present work provides 
the first test, for a 2d disordered particle system, of the 
full crossover between short- and long-distance behavior 
predicted [|l4| by conformal field theory. 

The values 0.54±0.01 and 0.57±0.05 we obtain for the 
exponent X t deviate somewhat from the value X t = 0.61 
found in previous work ^j. We attribute this difference 
mainly to the much larger statistical uncertainties in the 
numerical analysis of g. We also emphasize that the 
accuracies given here do not include systematic errors. 
Possible sources of these are the finite values of W/L 
and a/W, which affect the results obtained from different 
contact distances in a different way. We believe that the 
deviation between the values obtained from the moments 
(0.54) and the typical conductance (0.57) gives a realistic 




FIG. 2. The logarithm of the typical two-point conduc- 
tance (InT) of the Chalker-Coddington network at critical- 
ity versus the logarithm of the distance In |as | , for widths 
W = 3,4, 6, 9, 14, 20, and 30 (from bottom to top). The solid 
lines are functions given by Eq. ([H]). Inset: Rescaled curves 
(InT) + X t In W versus ln[x/W| with X t = 0.57. 

estimate of these systematic errors. 

R.K. thanks H. Moraal for help with the inversion of 
power series. This research was supported in part by the 
DFG, Sonderforschungsbereich 341. 



[1] E.R. Mucciolo, R.A. Jalabert, and J.-L. Pichard, J. Phys. 

I (France) 7, 1267 (1997). 
[2] E. Doron and U. Smilansky, Nonlinearity 5, 1055 (1992). 
[3] E.B. Bogomolny, Nonlinearity 5, 805 (1992). 
[4] J.T. Chalker and P.D. Coddington, J. Phys. C 21, 2665 

(1988). 

[5] M. Janssen, M. Metzler, and M.R. Zirnbauer, Phys. Rev. 

B 59, 15836 (1999). 
[6] Owing to the exchange symmetry T; m = T m i the l.h.s. of 

([]) can be symmetrized with respect to I and m, yielding 

■Kv(plf{pi/pm)+p m f{p m /pl)) c = (F(Ttm)). 

[7] M.R. Zirnbauer, Ann. Physik 3, 513 (1994). 
[8] J.T. Chalker and G.J. Daniell, Phys. Rev. Lett. 61, 593 
(1988). 

[9] This can always be done, as is shown by the following 
argument: On the closed network we choose an arbitrary 
closed path containing link m. Then we vary continu- 
ously the scattering parameters at the nodes lying on it 
so that the path becomes a simple loop insulated from the 
remaining network. A straightforward calculation shows 
that for this modified network adjusting ip m always suf- 
fices to tune one level to zero. By reasons of continuity 
and periodicity, this property persists as the nodes on the 
loop are put back to their initial values and the original 
network is restored. 
[10] R.E. Borland, Proc. Roy. Soc. London A 274, 529 (1963); 
A. Crisanti, G. Paladin, A. Vulpiani, Products of Random 
Matrices (Springer, Berlin, 1993). 



4 



[11] A.V. Kolesnikov and K.B. Efetov, Phys. Rev. Lett. 83, 
3689 (1999). 

[12] H. Schomerus and C.W.J. Beenakker, Phys. Rev. Lett. 

84, 3927 (2000). 
[13] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 
[14] J.L. Cardy, J. Phys. A 17, L385 (1984). 
[15] J.T. Chalker, J. Phys. C 21, L119 (1988). 
[16] A. Dohmen, P. Freche, and M. Janssen, Phys. Rev. Lett. 

76, 4207 (1996). 



5 



